Entropy production due to coupling to a heat 
bath in the kicked rotor problem 

Sankhasubhra Nag* Avijit Lahiri^and Gautam Ghosh* 
February 1, 2008 

Abstract 

Considering a kicked rotor coupled to a model heat bath both the clas- 
sical and quantum entropy productions are calculated exactly. Starting 
with an initial wave packet, the von Neuman entropy as a function of 
time is determined from the reduced density matrix while the Liouville 
evolution of the corresponding Husimi distribution provides us with the 
classical entropy. It is found that both these entropies agree reasonably 
satisfying the same asymptotic growth law and more importantly both 
are proportional to the classical Liapounov exponent. 

PACS numbcrs:03.65.Bz,05.40.+j,05.45.+b keywords: quantum open sys- 
tem, entropy 

The question of quantum classical correspondence in the context of open 
systems is an interesting one. It has been argued, for instance by Zurek and 
Paz[Q, that even a weak interaction with a random environment leads to the 
loss of quantum coherence and that classical behaviour may then be an emer- 
gent property of such systems. In other words, environmental effects may have 
a significant role to play in quantum classical correspondence in general and 
in the quantum dynamics of classically chaotic systems in particular. The von 
Neumann entropy production of an unstable oscillator interacting with a heat 
bath has been shown Q to depend monotonically on the classical local Lya- 
punov exponent. A more extensive study of the kicked rotor in a bath by Miller 
and Sarkar || brings out in greater detail this correspondence between quantum 
entropy production and classical Lyapunov exponent. The quantum evolution 
of classically chaotic open systems has been shown to possess a sensitive depen- 
dence on initial conditions || and provides another interesting indication of the 
role of the environment. As a further instance one may mention the high degree 
of sensitivity of quantum dynamic localization Q to the presence of noise found 
in a study on the kicked rotor ||. Extensive studies of the kicked rotor coupled 
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to a bath were carried out by Dittrich and Graham^ and by Cohen|Q. These 
authors were interested in reducing the problem and showing its equivalence to 
a classical stochastic map in the semiclassical limit with the aid of the Wigner 
function representation of the density matrix. 

In the present paper we look into a few aspects of quantum classical cor- 
respondence in open systems by focussing on the kicked rotor interacting with 
a bath of harmonic oscillators through what is termed a nondemolition cou- 
pling. The latter involves an interaction Hamiltonian commuting with the sys- 
tem Hamiltonian as a result of which the reduced density matrix can be obtained 
exactly This feature has been used in in computing the von Neumann 
entropy production rate of the kicked rotor. We show that the special nature 
of the coupling allows one to obtain the classical reduced distribution function 
[^) as well and then to compare the diffusion rate and the entropy production 
at the classical and quantum levels under similar initial conditions. Numerical 
computations show that in the semiclassical limit the asymptotic von Neuman 
entropy follows the same classical growth law as obtains for a diffusion process 
viz. s(t) = A + B\nt and that this regime is established smoothly after a very 
few kicks even for relatively weak bath couplings. Of more interest is the fact 
that in the asymptotic region the quantum entropy is proportional to the clas- 
sical Liapounov exponent which may perhaps be looked upon as a signature of 
quantum chaos in semiclassical quantum mechanics. 

We take a cylindrical phase space (0 < q < 2ir; — oo < p < oo) and define 
the kicked rotor Hamiltonian to be, 

+00 

Hrotor = p 2 /2 + Kcosq 5{t-n). (1) 

n— — oo 

The classical dynamics is defined by the standard map, 

Pn+i = Pn + Ksmq n+1 , (2) 
q n +i = q n + Pn mod 2ir, (3) 

where q n and p n are the values just after the nth kick. Below a critical kick 
strength K c (^ .97) KAM tori prevent the indefinite growth in energy while 
above K c , one has a mixed phase space of islands surrounded by a chaotic sea 
and the diffusive growth law holds viz. {{p n —Po) 2 ) = D(K)n, where the average 
is over a distribution of initial points and the diffusion coefficient tends to K 2 /2 
for large K. In quantum mechanics the evolution is governed by the single step 
unitary evolution operator which takes the state from just after the nth kick to 
just after the (n + l)th kick and is given by, 

iK 

U(n + l,n) = U k U f = ex P[-^- cos l\ exp[-ip 2 / (2h)}. (4) 

In the momentum representation Uf is just a phase factor and the matrix 
element for Uk is given by, 

(l\U k \m) = (i) l - m J^ m (K/h). (5) 
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Numerical computation with an initial wave packet shows a diffusive growth in 
energy for a finite time determined by the values of K and h after which the 
energy oscillates quasiperiodically about an average value ||, [h| . 

The total Hamiltonian for a system with a nondemolition coupling to an 
oscillator bath is chosen to be, 

k 

+<f>(H ) Y^c k ^T k q k + l - V^-0 2 (iJ o ), (6) 
fe 2 k " k 

where the first two terms are the system and bath Hamiltonians respectively. 
The third term is the interaction with c k being the coupling constant to the fcth 
mode, 4> being an arbitrarily chosen function, while the fourth one is a renor- 
malization term. We will work out the time evolution for a general Hamiltonian 
H (p) and then specialize to the case of the rotor. Notice that the interac- 
tion is a function of the system Hamiltonian and therefore momentum will be 
conserved. Making a canonical transformation with the generating function, 

F(q, q k ; P, P k ) = qP + Y, IkPk + 4>(P) £ -^jflk, (7) 
fc fc w fc 

the Hamiltonian is reduced to the uncoupled form, 

H = H {P) + \^ P k+^lQl)- (8) 

fc 

Solving the equations of motion obtained from (||) and reverting back to orig- 
inal variables we get our required solutions viz. q(0) and p(0) in terms of 
q,p,qk(0) andpfc(O). As an initial point {<?(()), p(0)} evolves in time t to {q,p} 
the distribution function evolves from f°[q(0),p(0)] to f[q,p] maintaining its 
value at the corresponding points. However, since the evolution depends on 
the bath variables, the reduced distribution function is obtained by integrating 
over the bath quantities which are assumed initially to have values pertaining 
to thermal equilibrium. Thus, 

fred.{q,P,t) = f f[q(0),p(0)]]l—dpk(0)dq k (0), (9) 

J fc Zk 

where in the r.h.s. [g(0),p(0)] are understood to be written in terms of [q,p\. 
The factor g k is the thermal probability for the fcth bath mode at temperature 
T(= 1/p) and it is normalized by the partition function z k . We carry out the 
integration in equation (|^) and go to the continuum limit whereby the sum 
over modes is replaced by an integration over frequency with an ohmic spectral 
density function given by, 

£ - J f(u)du, /M = v ^-je-"/«°. (10) 
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The reduced distribution function is given by, 



M*,P,t) = y^ydrn ff -(-^)* + r 

/ fir 2 
-</)0'?7arctanw c t,p] exp(-^-^y), (11) 

s(t) = 2r/0' 2 [iarctanw c i- - — ln(l + w^t 2 )]. (12) 

2u c 

If we look at the fundamental solution of ( |ll| ) when /° is a 5-function we find 
/red. to be a gaussian . For uo c t <C 1, s(i) oc i 2 and the gaussian spreads fast 
but after some time depending on u> c , arctanw c i saturates at ir/2 and s(t) oc t, 
so that equation ( pi] ) describes a normal diffusion process. In the case of the 
rotor (dH n (p)/dp) = p, and for the coupling we take <f)(p) — p, so that (j)' = 1. 
The fourth term in (^) therefore renormalizes the rotor mass. Because of the 
nondcmolition nature of the coupling the p-distribution remains the same in 
the absence of the kicking term whereas the g-distribution suffers a shift and 
a gaussian local averaging. However, as it turns out, the bath brings about a 
mixing in the p-distribution on a fine scale when the kicks couple q and p. . 

For the quantized form of the Hamiltonian (^|) the matrix elements of the 
reduced density operator for the rotor at time t are given by |^|. || , 

(m | p re d.{t) I n) — exp[— ih(m 2 — n 2 )t/2] 

(m | p(0) | n) eM-i^n ~ A 2 mn ], (13) 



where, 



A mn = VH™ 2 ~ n 2 ) tan" 1 w c t, (14) 
A 2 mn = ^(m-n) 2 ln(l+^ 2 i 2 ) 

oo 
k—1 



Notice that the diagonal matrix elements are unaffected by the interaction with 
the bath when the kicking term is absent. 

In analysing the quantum evolution we take an initial normalized wave 
packet of the form, | ip) — a; | I) which is localized and peaked around 
the momentum p = HIq- We achieve similar initial conditions for the classical 
evolution by taking the Husimi distribution JTT| corresponding to | ip). This is 
given by, 



f(q,p)=\{q,p\i>)\ 2 = J2^peM-P 2 /f> 2 



rn . n 



+ (m + n)p/h - (m + n )/2] cos(m - n)q, (16) 



4 



which is a diagonal approximation of the density matrix in the coherent state 
representation. f(q,p) is periodic in q with period 2tt so that ( [l6| ) is really 
defined on the cylinder. The coefficients ai are chosen such that | tp) represents 
a wave packet with Ap = Aq = h/2. This is achieved by parametrising a; in 
the form, a; = N exp[— al 2 + bl], where a is real and b is complex such that the 
three parameters a and the real and imaginary parts of b determine the centre 
and the width of the wave packet. 

For the quantum evolution the computation starts with the density matrix 
corresponding to | ip) . After unit time step the reduced density matrix is evalu- 
ated using equation ([l3|). The unitary kick operator Uk then connects the matrix 
elements for the density operators just before and after the kick according to 
(4). The von Neumann entropy is calculated by using the definition, 

S = —Trp In p = — Xj In Xj , (17) 

i 

where the Ai are the eigenvalues of the density operator. The computation is 
repeated till the entropy reaches its asymptotic growth rate. To compute the 
classical evolution we start with the Husimi distribution corresponding to the 
initial quantum wave packet and use equation ([)]) to find the reduced distribu- 
tion after unit time step. The distributions just before and after the kicks are 
related by a simple shift in argument viz. f^ ed (q,p) — f~ ed (q,p — Ksiaq). The 
entropy is calculated by using the formula, 

./ 2ttH I 2nh y ' 



It is known that Sh calculated from f(q,p) as defined in (16) satisfies the in- 
equality, Sh > S. It turns out that the classical approximation is good as long 
as f(q,p) is a smooth function spread over an area in phase space that is much 
greater than 2irh. If there are small distance fluctuations or if f{q,p) is concen- 
trated on a small region of the phase space, then the classical approximation 
can be very bad |J. This fact will be borne out in our numerical computations. 
Notice that the von Neumann entropy for the initial pure state is zero. However, 
this is not true of (|is|). The Husimi distribution corresponding to the initial 
wave packet is sufficiently well localised in q to be considered a gaussian with 
ApAq — h and its entropy is just unity.The distribution is chosen to be centred 
at q = 0,p = irH which lies in a chaotic region for all the K- values cosidered by 
us. 

To understand the diffusion process in detail we compute the (/-distribution, 
gi(q) and the p-distribution, gi (p) as the integrals of f(q,p) over p and q re- 
spectively. The diffusion rate in q is controlled by the parameter lo c and (3. For 
our choice, uj c — 5/S« 10 and (3 = 0.1 satisfying [3ftw c < 1, we find that gi(q) 
almost reaches the value (^) before the first kick [fig. (|]a)]. Decreasing w c , will 
increase the equilibriation time for the g-distribution. The large scale structure 
of the p-distribution on the other hand is determined by the kicks, the bath 
providing a mixing on a finer scale through the coupling with the kicks. Figure 
(0b) shows the initial evolution of #2 (p) ■ 
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Figure 1: Evolution of (a) g\{q) and (b) ^(f) f° r V — 1-0, ^ = 3.5. 
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Figure 2: Comparison of classical and quantum energy growth for K = 3.5, H 
0.46. 



Next, we compute the classical and quantum energies. The effect of h for 
fixed K (the relevant parameter actually is K/K) on the variation of (E) with t 
as also on the entropy production (see below; neither of these two is shown in our 
figures) indicates that one can distinguish between a semiclassical and strongly 
quantum regime: while oscillations due to quantum correlations are pronounced 
in the latter, they are almost absent in the former. All our numerical work 
relates mainly to the semiclassical regime with H taken as a rational multiple 
of the golden mean (v5 — l)/2 to avoid resonances. The energy growth is 
typically diffusive for both classical and quantum dynamics and the diffusion 
rate increases with K. What is interesting is the effect of the bath coupling 
strength rj on the energy growth (fig.^): on increasing rj from a low value one 
finds that the classical and quantum energy growth curves quickly approach 
each other (cf. Ott, Antonsen and Hanson and their growth rates saturate 
at K 2 /4 even for comparatively low values of K. This is consistent with the 
role of the bath outlined above. 

Figures (||(a,b)) present results on entropy production. Note that there is 
an h dependence in the classical entropy Sh arising from the Husimi density 
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Figure 3: Entropy growth in the quantum case for (a) different h with K = 
3.5, r\ = 1(T 2 and (b) different r) with K = 3.5, H = 0.46. 



corresponding to an initial wave packet. The evolution of the Husimi density 
satisfies the Liuville equation in the lowest order in h and that is what we have 
considered here. The classical and quantum entropy productions are compared 
in fig. |3|(a) which shows that the two converges quickly excepting for compara- 
tively large values of H and even when they differ their asymptotic growth rates 
agree. The time dependence of entropy production is elucidated in fig. ||(b) 
which uses a logarithmic time scale: asymptotically the entropy settles down to 
a logarithmic growth. Curves for different values of rj once again show a quick 
saturation with increasing rj. 

The asymptotic entropy production fits nicely with the formula S = A + 
Bhin which we now explain by referring to the classical entropy Sh- Entropy 
production occurs due to the coarse graining provided by the bath along with the 
diffusion provided by the kicks. The bath quickly uniformizes the q-distribution 
and converts the p-distribution to a smoothed gaussian so that in the asymptotic 
regime 

V27rAp 2Ap^ 

Together with the diffusive growth law, Ap — {K/y/2)n^ this yields for the 
entropy, 

S H = ^+H^-)+\nK+hnn. (19) 

Since for large K, the Lyapunov exponent is, A = hi-^, the coefficient A is 
linear in the Lyapunov exponent for large K. 

The growth law S = A + B In n with A depending linearly on In K and 
B = 1/2 is numerically found to apply to the quantum entropy as well in the 
semiclassical regime. Figures ||(a,b) depicts the variation of A with In if as 
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Figure 4: 
quantum 
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Variation of A with In K for different r\ in the (a) classical and (b) 
case, (c) shows a comparison with the values estimated from (|l^) for 
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Plot of B with K for different ?/ in (a)classical case and (b)quantum 



computed numerically from the classical and quantum entropy production data 
where the linear dependence is evident and where the role of r\ in quantum 
classical correspondence is apparent once again. Figure ^(c) compares the nu- 
merically computed values of A with the estimated values from equation (19). 

Figures |^(a,b) present corresponding data for the coefficient B which is in- 
deed found to be close to 1/2. The discrepancy is explained by the presence 
of residual phase space barriers to diffusion as also incomplete mixing which 
disappears more and more with increasing K and rj. In figures ^ and h is 
held at 0.46. 

In conclusion, the findings presented above underline the role of the bath 
in establishing a quantum-classical correspondence wherein the quantum en- 
tropy production carries with it the characteristics of the classical phase space 
structure: for instance, the quantum entropy is tied to the classical Lyapunov 
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exponent |y, g|. It would be interesting to look into the entropy production in 
the kicked Harper model where the phase space is compact and the diffusive 
growth law (Ap oc y/n) does not hold. 
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